Epitaxial strain effects in the spinel ferrites CoFe204 and NiFe204 from first principles 
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The inverse spinels CoFe2 04 and NiFe2 04, which have been of particular interest over the past 
few years as building blocks of artificial multiferroic heterostructures and as possible spin-filter 
materials, are investigated by means of density functional theory calculations. We address the 
effect of epitaxial strain on the magneto-crystalline anisotropy and show that, in agreement with 
experimental observations, tensile strain favors perpendicular anisotropy, whereas compressive strain 
favors in-plane orientation of the magnetization. Our calculated magnetostriction constants Aioo of 
about —220 ppm for CoFe204 and —45 ppm for NiFe204 agree well with available experimental 
data. We analyze the effect of different cation arrangements used to represent the inverse spinel 
structure and show that both LSDA+J7 and GGA-I-J7 allow for a good quantitative description of 
these materials. Our results open the way for further computational investigations of spinel ferrites. 

PACS numbers: 



I. INTRODUCTION 

Spinel ferrites CoFe204 and NiFe204 are insulating 
magnetic oxides with high magnetic ordering temper- 
atures and large saturation magnetizations. -'^ This rare 
combination of properties makes them very attractive 
for a wide range of applications. Recently, particular 
attention has been focused on the possible use of spinel 
ferrites as magnetic components in artificial multiferroic 
heterostructures^^— or as spin-filtering tunnel barriers for 
spintronics devices i^^— 

For these applications, the corresponding materials 
have to be prepared either in the form of thin films, 
grown on different substrates, or as components of more 
complex epitaxial heterostructures .i^"— Due to the mis- 
match in lattice constants and thermal expansion coeffi- 
cients between the thin film material and the substrate, 
significant amounts of strain can be incorporated in such 
epitaxial thin film structures, depending on the specific 
growth conditions and substrate materials. This epitax- 
ial strain can then lead to drastic changes in the proper- 
ties of the thin film material. Indeed, a reoriention of the 
magnetic easy axis under different conditions has been 



reported for CoFe204 



and a strong enhancement 



of magnetization and conductivity has been observed in 
NiFe2 04 thin films J^^i^i^ 

In order to efficiently optimize the properties of thin 
film materials, it is important to clarify whether the ob- 
served deviations from bulk behavior are indeed due to 
the epitaxial strain or whether they are induced by other 
factors, such as for example defects, off-stoichiometry, 
or genuine interface effects. First principles calculations 
based on density functional theory (DFT),— can pro- 
vide valuable insights in this respect by allowing to ad- 
dress each of these effects separately. 

Here we present results of DFT calculations for the 
structural and magnetic properties of epitaxially strained 
CoFe204 and NiFe204, with special emphasis on strain- 
induced changes in the magneto-crystalline anisotropy 
energy (MAE). Our results are representative for (001)- 



oriented thin films of CoFe204 and NiFe2 04, grown on 
different lattice-mismatched substrates. Our results pro- 
vide important reference data for the interpretation of 
experimental observations in spinel ferrite thin films and 
in heterostructures consisting of combinations of spinel 
ferrites with other materials, such as perovskite struc- 
ture oxides. 

We find a large and strongly strain-dependent MAE 
for CoFe204, and a smaller but also strongly strain- 
dependent MAE for NiFe2 04. We discuss the infiuence 
of different cation arrangements within the inverse spinel 
structure and analyze the difference in the structural and 
magnetic properties due to different exchange-correlation 
functionals used in the calculations. From our calcula- 
tions we obtain the magnetostriction constants Aioo for 
both CoFe204 and NiFe204, which agree well with avail- 
able experimental data. 

This paper is organized as follows. Sec. Ill Al gives 
a brief overview over the properties of CoFe204 and 
NiFe2 04 that are important for the present work and also 
summarizes results of previous DFT calculations. The 
basic equations governing the magnetoelastic properties 
of cubic crystals are presented in Sec. IIIBl followed by 
a detailed description of the structural relaxations per- 
formed in this work in Sec. Ill C[ and a summary of fur- 
ther computational details in Sec. Ill D[ The results of the 
bulk structural properties will be presented in Sec. IIII Al 
whereas the effect of strain on the structural properties 
is analyzed in Sec. IIIIBI The effect of strain on the MAE 
is discussed in Sec. IIII CI Finally, in Sec. IIVI a summary 
of our main conclusions is given. 



II. BACKGROUND AND COMPUTATIONAL 
DETAILS 

A. Spinel structure and previous work on ferrites 

The spinel structure (space group Fd3m, general for- 
mula AB2X4) contains two inequivalent cation sites, the 
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tetrahedrally-coordinated A site {Td symmetry, Wyck- 
off position 8a), and the octahedrally-coordinated B site 
(Oh symmetry, Wyckoff position 16d). In the normal 
spinel structure, ah A sites are occupied by one cation 
species (divalent cation), whereas ah B sites are occu- 
pied by the other cation species (trivalent cation). On 
the other hand, in the inverse spinel structure the triva- 
lent cations occupy all A sites as well as 50 % of the 
B sites, whereas the remaining 50 % of the B sites are 
occupied by the divalent cations. If the distribution of di- 
valent and trivalent cations on the B sites is completely 
random, all B sites remain crystallographically equiva- 
lent and the overall cubic FdZm symmetry is preserved 
(see Fig. 1(a) ). 

Both CoFe204 and NiFe2 04 crystallize in the inverse 
spinel structure, even though for CoFe204 the inversion is 
typically not fully complete, i.e. there is a nonzero Co^+ 
occupation on the A site. Thereby, the exact degree of in- 




(a) 




(b) 



FIG. 1: (Color online) Inverse spinel structure of spinel fer- 
rites. Oxygen cations are depicted as small red spheres, 
and the coordination polyhedra surrounding cation sites are 
shaded, (a) Random cation distribution on the octahedrally 
coordinated B sites, i.e. all B sites remain equivalent, (b) 
Imma cation arrangement used throughout this work. In- 
equivalent B sites are represented by different shadings of the 
corresponding octahedra. 



version depends strongly on the preparation conditionsi^ 
To the best of our knowledge, no deviations from cu- 
bic symmetry have been reported for either system, i.e. 
the Co^"'"/Ni^+ cations are believed to be randomly dis- 
tributed on the B sites/22, 

According to the formal and electron configura- 
tions corresponding to Fe^+ and Co^+, respectively, these 
ions can in principle exhibit both high-spin and low-spin 
states, but only the high-spin states are experimentally 
observed in both spinel ferrites. For the electron con- 
figuration of Ni^+ no such distinction exists. The mag- 
netic moments of the A site cations are oriented antipar- 
allel to the magnetic moments of the B site cations: the 
so-called Neel-type ferrimagnetic arrangement Ji Thus, 
the magnetic moments of the Fe'^+ cations on the A and 
B sites cancel each other exactly, and the net magneti- 
zation is mainly due to the divalent B site cations, i.e. 
either Co^"*" or Ni^"*". This results in a magnetic moment 
per formula unit close to the formal values of 3 /.ts and 
2 /is in CoFe204 and NiFe204, respectively. Deviations 
from these values can be either due to orbital contribu- 
tions to the magnetic moments, or due to incomplete 
inversion and off-stoichiometric cation distribution. 

Both CoFe204 and NiFe204 are small gap insulators, 
but information on the experimental gap size is very lim- 
ited. Waldron used infrared spectra to obtain threshold 
values of 0.11 eV and 0.33 eV for the electronic transi- 
tions in CoFe204 and NiFe204, respectively?^ whereas 
Jonker estimated the energy gap in CoFe204 to be 
0.55 eV, based on resistivity measurements along with 
other methods.^ 

Theoretical calculations of the electronic structure of 
spinel ferrites so far have been focused mostly on mag- 
netite (Fe304). This material can be viewed as parent 
compound for the spinel ferrites, including CoFe204 and 
NiFe2 04, which are obtained by substituting the Fe^"*" 
cation in magnetite by a different divalent 3(i transition 
metal cation. In an early work, Penicaud et al. per- 
formed DFT calculations within the local spin-density 
approximation (LSDA) for magnetite and the respective 
Co-, Ni-, Mn-, and Zn-substituted ferritesj24 The use of 
LSDA leads to half-metallic band-structures for all sys- 
tems except NiFe2 04, in contrast to the insulating char- 
acter observed experimentally. (We note that the case 
of magnetite is somewhat more involved than that of the 
other spinel ferrites, since magnetite exhibits a metal- 
insulator transition at ~120 K.) It was later shown by 
Antonov et al. that insulating solutions for Co-, Ni-, 
and Mn- substituted Fe3 04 can be obtained within the 
LSDA-I-C/ approach. The same was found by Szotek et 
al. using a self-interaction-corrected LSDA approachi^ 
The latter study also addressed the energetic difference 
between normal and inverse spinel structures with dif- 
ferent valence configurations. The electronic structure 
of NiFe204 was also calculated within a hybrid func- 
tional approach, where a large band-gap of 4 eV was 
obtained by using 40 % of Hartree-Fock exchange in the 
exchange-correlation energy functional. Recently, Per- 



3 



ron et al. investigated different magnetic arrangements 
for NiFe204 in both normal and inverse spinel structures 
using both LSDA and the generalized gradient approx- 
imation (GGA), and found the inverse spinel structure 
with Neel-type ferrimagnetic order to be energetically 
most favorable,— in agreement with the experimental ob- 
servations. 

Calculations of the MAE in strained CoFe204 and 
NiFe2 04 have been reported by Jeng and Guoi^i^ How- 
ever, due to the use of the LSDA, these calculations were 
based on half-metallic band-structures for both materi- 
als. Furthermore, no information on structural properties 
or the influence of the specific cation arrangement used 
in the calculation were given. 



B. Magnetoelastic energy of a cubic crystal 

In this work we are concerned with the effect of epi- 
taxial strain on the structural and magnetic properties 
of CoFe204 and NiFe204, i.e. with the elastic and mag- 
netoelastic response of these systems. Here, we therefore 
give a brief overview over the general magnetoelastic the- 
ory for a cubic crystal, and present the most important 
equations that are used in Sec. IHII to analyze the results 
of our first principles calculations. 

The magnetoelastic energy density / — E/V of a cubic 
crystal can be written asi^ 



f — Ik + fcl + /mcl : 



(1) 



where the three individual terms describe the cubic (un- 
strained) magnetic anisotropy energy density /k, the 
purely elastic energy density /d and the coupled magne- 
toelastic contribution /moi, respectively. To lowest order 
in the strain tensor £y and in the direction cosines of 
the magnetization vector, these terms have the following 
forms: 
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Here, K denotes the lowest order cubic anisotropy con- 
stant, Cii, Ci2, and C44 are the elastic moduli, and Bi 
and B2 are magnetoelastic coupling constants. 
The bulk modulus B is defined as: 



B^Vo 



(5) 



where i^tot is the total energy and Vq is the equilibrium 
bulk volume. Using Eqs. ((l|)-(|4]) the bulk modulus of 
a cubic crystal can be expressed in terms of the elastic 
moduli Cii and C12: 



S = i {Cn + 2C12) . 



(6) 



In this work we investigate the effect of epitaxial strain 
that is induced in thin film samples by the lattice mis- 
match to the substrate. This situation can be described 
by a fixed in-plane strain Exx = £yy = (ct ~ olq)/o,Oi where 
a is the in-plane lattice constant of the thin film material 
and Oq is the corresponding lattice constant in the bulk. 
The resulting out-of-plane strain e^z = (c — ao)/ao can be 
obtained from Eqs. ©-(Ill) together with the condition of 
vanishing stress for the out-of-plane lattice constant c, i.e. 

= 0. In the demagnetized state e^z is related to the 
applied in-plane strain via the so-called two-dimensional 
Poisson ratio V2d'-^^ 



V2D 



^ ZZ 



(7) 



After calculating both z^2D and the bulk modulus using 
DFT, the elastic moduli Cn and C12 can thus be ob- 
tained from Eqs. (jG]) and Q. 

In Sec. nil CI we monitor the differences in total en- 
ergy for different orientations of the magnetization as a 
function of the in-plane constraint Exx- Using expres- 
sion ^ for 1^20 and taking the energy for orientation of 
the magnetization along the [001] direction as reference, 
i.e. Afhu = /ooi - fhkh one obtains: 

^/(lOO) = -Bi{l'2D + 1) £xx , 
A/(110) = -Bi{v2D + 1) txx - , 

2 3 (8) 

= -^Bi{v2D + 1) ^xx - -^K , 

A/(ioi) = -]^Bi{v2D + l)exx - . 

Thus, the epitaxial strain-dependence of these energy dif- 
ferences is governed by the magnetoelastic coupling con- 
stant Bi and the two-dimensional Poisson ratio 1^20 ■ 

Since the constant Bi is not directly accessible by 
experiment, the linear magnetoelastic response is typ- 
ically characterized by the magnetostriction constant 
Aioo, which is related to Bi and the elastic moduli Cn 
and C12: 



Bi 



3 Cn — C12 



(9) 



{V=Vo) 



Aioo characterizes the relative change in length (lattice 
constant) along [100] when the material is magnetized 
along this direction, compared to the unmagnetized state. 



structural relaxations for the inverse spinel 
structure 



As described in Sec. Ill Al the distribution of divalent 
and trivalent cations on the octahedrahy coordinated B 
sites in the inverse spinel structure is assumed to be 
random for both CoFe2 04 and NiFe2 04. On the other 
hand, the periodic boundary conditions employed in our 
calculations always correspond to a specific cation ar- 
rangement with perfect long-range order. Even though 
a "quasi-random" distribution of divalent and trivalent 
cations could in principle be achieved by using a very 
large unit cell, the required computational effort would 
be prohibitively large. For simplicity, we therefore re- 
strict ourselves to the smallest possible unit cell which 
contains two spinel formula units, i.e. four B sites. Dis- 
tributing two Fe atoms on two of these sites, and filling 
the other two sites with either Co or Ni, lowers the sym- 
metry from spacegroup FdZm (#227) to Imma (#74) 
independent of which of the four B sites are occupied by 



Fe (see Fig. 1(b)). Different choices ("settings") simply 



lead to different orientations of the orthorhombic axes 
relative to the Cartesian directions. It will become clear 
from the results presented in Sec. IIIII that the specific 
cation arrangement used in our calculations does not crit- 
ically affect our conclusions. 

Within the lower Imma symmetry (and setting 1, see 
below), the tetrahedrally coordinated A sites are located 
on Wyckoff position 4e (0, |,z), whereas the octahe- 
drahy coordinated B sites split into Wyckoff positions 
4b (0, 0, ^) and 4d (i, j, |). In addition, the oxygen po- 
sitions split into Wyckoff positions 8i {x,j,z) and 8h 
(0,y,^). 

In order to minimize the effect of this artificial sym- 
metry lowering, and to obtain results that are as close 
as possible to the average cubic bulk symmetry seen in 
experiments, we apply the following constraints during 
our structural relaxations. For the calculations corre- 
sponding to the unstrained bulk case, we constrain the 
lattice parameters along the three cartesian direction to 
be equal, a — b = c, and we fix the A site cations to their 
ideal cubic positions, corresponding to 2;(4e) = |. On 
the other hand, since the oxygen positions are charac- 
terized by one free structural parameter already within 
cubic Fd3m symmetry (Wyckoff position 32e {u,u,u)), 
we do not apply any constraints to the 8i and 8h posi- 
tions within Imma symmetry. The relaxed bulk struc- 
ture is then found by relaxing all internal positions for 
different volumes and finding the volume that minimizes 
the total energy. For the relaxations corresponding to 
a certain value of epitaxial strain, we constrain the two 
in-plane lattice constants to be equal, and then vary the 
out-of-plane lattice constant c while relaxing all internal 
coordinates (except for z(4e) = |). 

While for the unstrained bulk structure all possible 
B site cation distributions that can be accommodated 
within the primitive fee unit cell of the spinel struc- 
ture lead to the same Imma symmetry, different cases 



z=o 



z=0 




FIG. 2: Cation distribution within different x-y planes for dif- 
ferent value s of z, corresponding to "setting" 1 (a) and "set- 
ting" 3 (b) (see main text for details). Only layers contain- 
ing B site cations are shown. Divalent cations (Co^''', Ni^"'") 
are depicted as filled black circles, trivalent cations (Fe^^) as 
open black circles, and oxygen anions as black crosses. Note 
that we use the convention where the origin is located at the 
midpoint between two A sites. 



can be distinguished once an epitaxial constraint is ap- 
plied to the resulting structure. To illustrate this, the 
cation distributions for two different settings are depicted 
in Figs. 2(a) and |2(b)[ The B site cation distribution 
(and the oxygen positions) within different x-y planes 
are shown, corresponding to different "heights" z. The 
different settings simply correspond to different orienta- 
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tions of the orthorhombic Imma crystal axes relative to 
the original cubic axes. 

It can be seen from Fig.[2]that the B sites are arranged 
in an interconnected network of chains along (llO)-type 
directions. In setting 1 (Fig. 2(a) I, the corresponding 
chains within the same x-y plane contain alternating di- 
valent and trivalent cations. In contrast, for setting 3 



(Fig. 2(b) I each x-y plane contains only one unique cation 



species, which then alternates between adjacent planes 



along the z direction (see also Fig. 1(b)). In setting 1 
the same alternating planes are oriented perpendicular to 
the X direction, whereas for setting 2 (not shown) these 
planes are oriented perpendicular to the y direction. If an 
epitaxial constraint is applied within the x-y plane, set- 
ting 3 becomes different from settings 1 and 2, with the 
difference being the orientation of the "substrate plane" 
relative to the planes defined by the cation order. 

In view of this, we have performed all calculations cor- 
responding to epitaxially strained systems for both set- 
ting 1 and setting 3. The differences between the results 
obtained for the two different settings then represent a 
measure for the sensitivity of these results from the spe- 
cific cation arrangement used in our calculations. 



D. Other computational details 



TABLE L Optimized bulk lattice constants ao and bulk 
moduli B for CoFe2 04 and NiFe2 04, calculated using the 
LSDA-|-)7, GGA, and GGA-|~f7 exchange-correlation func- 
tionals in comparison with experimental data. 





CoFe 


2O4 


NiFe2 04 




ao [A] 


B [GPa] 


ao [A] 


B [GPa] 


LSDA+(7 


8.231 


206.0 


8.196 


213.1 


GGA 


8.366 


211.0 


8.346 


166.2 


GGA+;7 


8.463 


172.3 


8.426 


177.1 


Exp. (Ref._42) 


8.392 


185.7 


8.339 


198.2 


Exp. (Ref. 43) 


8.35 




8.325 





UcS — U ~ J — 3 eV is used for all transition metal 
cations. The corresponding results are compared to pure 
GGA calculations, using the GGA approach of Perdew, 
Burke, and Ernzerhof4i (We restricted the comparison to 
pure GGA since Perron et al^ presented some evidence 
(for NiFe204) that LSD A might not be appropriate to 
properly describe these materials.) 

Values for the local magnetic moments and atom- 
projected DOS are obtained by integration of the ap- 
propriate quantities over atom-centred spheres with radii 
taken from the applied PAW potentials (1.164 A (Fe), 
1.302 A (Co), and 1.058 A (Ni)), respectively. 



All calculations presented in this work were performed 
using the projector-augmented wave (PAW) method^^ 
implemented in the Vienna ab initio simulation pack- 
age (VASP 4.6),^""^ Standard PAW potentials supplied 
with VASP were used in the calculations, contributing 
nine valence electrons per Co (4s^3d^), 16 valence elec- 
trons per Ni (3p^4s^3d^), 14 valence electrons per Fe 
(3p<'4s23d^), and 6 valence electrons per O {2s'^2p^). A 
plane wave energy cutoff of 500 eV was used, and the 
Brillouin zone was sampled using different fc-point grids 
centered at the F point. A 5 x 5 x 5 fc-point grid was 
used for the structural optimization and all total energy 
calculations, whereas a finer 7x7x7 grid was used 
to calculate densities of states (DOS). The tetrahedron 
method with Blochl corrections was used for Brillouin 
zone integratioui^ We have verified that all quantities of 
interest, in particular the MAEs, are well converged for 
the used fc-point grid and energy cutoff. All structural re- 
laxations were performed within a scalar-relativistic ap- 
proximation, whereas spin-orbit coupling was included in 
the calculation of the MAEs. 

As already noted in Sec. Ill Al CoFe204 and NiFe204 
are small gap insulators, whereas half-metallic or, 
in the case of NiFe2 04, only marginally insulating 
band-structures have been obtained in previous LSDA 
calculations.— In the present work we therefore use the 
LSDA-|-?7 and GGA-\-U approachj^^ which is known to 
give a good description of the electronic structure for 
many transition metal oxidesj^ We employ the Hubbard 
"-I-C/" correction in the simplified, rotationally invari- 
ant version of Dudarev et al.^ where the same value 



III. RESULTS AND DISCUSSION 
A. Unstrained bulk structures 

We first present our results for the unstrained bulk 
structures. The calculated lattice constants and bulk 
moduli for both CoFe204 and NiFe2 04 using different 
exchange-correlation functionals are summarized in Ta- 
bic H] It can be seen that for both CoFe204 and NiFe204 
the use of LSDA-|-[/ leads to an underestimation of the 
lattice constant and an overestimation of the bulk mod- 
ulus compared to the experimental values, whereas the 
opposite is the case for GGA-\-U . The corresponding de- 
viations (~l-2 % for the lattice constants) are typical for 
complex transition metal oxides (see e.g. Refs. I44l - l46l ). 

Interestingly, the lattice constants calculated within 
pure GGA match the experimental values almost per- 
fectly for both CoFe204 and NiFe2 04. However, this is 
somewhat fortuitous and probably due to a cancellation 
of errors, as can be seen by the large discrepancies in 
the bulk moduli. It will become clear in the following, 
that the "-I-C/" correction is necessary in order to obtain 
a good description of the electronic structure for both 
CoFe204 and NiFe204. 

As discussed in Sec. Ill Ci the cation arrangement in our 
unit cell lowers the symmetry to orthorhombic Imma, 
with 4 independent parameters describing the positions 
of the oxygen anions at Wyckoff positions 8h and Si, com- 
pared to one parameter for Wyckoff position 32e in the 
cubic space group Fdim. Table |TT] lists the correspond- 



TABLE II: Calculated Wyckoff parameters (setting 1) for 
the oxygen anions 81i {x,j,z) and 8i {0,y,z) for CoFe204 
and NiFe2 04 using different exchange-correlation function- 
als. The last line lists the corresponding parameters resulting 
from Wyckoff position 32e (u, u, u) within Fd3m symmetry, u 
is obtained from these relations by averaging over recalculated 
u values for each dataset. 



CoFe2 04 




8i 






8h 








X 




z 


y 




z 


u 


LDA+17 


0.235 




-0.498 


0.009 




-0.257 


0.255 


GGA 


0.240 




-0.496 


0.008 




-0.255 


0.255 


GGA+U 


0.234 




-0.499 


0.007 




-0.259 


0.255 


NiFe204 




Si 






8h 








X 




z 


y 




z 


u 


LDA+C/ 


0.237 




-0.495 


0.010 




-0.258 


0.256 


GGA 


0.239 




-0.496 


0.009 




-0.257 


0.255 


GGA+;7 


0.235 




-0.496 


0.008 




-0.258 


0.256 


Fdim 


3 






2^-1 




—u 





ing WyckofF parameters obtained from our structural op- 
timizations. It can be seen that differences between dif- 
ferent exchange-correlation functionals are rather small. 

The last line in Table [ll] indicates the relation be- 
tween the Wyckoff positions 8i and 8h in the Imma space 
group (setting 1) and Wyckoff position 32e {u, u, u) in cu- 
bic FdZm symmetry (assuming that no actual symmetry 
breaking occurs). These relations allow us to obtain an 
average Wyckoff parameter u, by calculating the value of 
u corresponding to each of the four calculated Wyckoff 
parameters a;(8i), z{%i), 2/(8h), and z(8h) for each data 
set and subsequent averaging. The resulting values for u 
agree very well with available experimental data, which 
are 0.256 for CoFe204 and 0.257 for NiFe204.— Further- 
more, the values calculated from the individual 8h and 
8i Wyckoff parameters deviate only very little from the 
average values, which indicates that the lower symmetry 
used in our calculation has only a negligible effect on the 
internal structural parameters. 

Figs. 12] and 0] show the calculated spin-decomposed 
densities of states (DOS) for CoFe204 and NiFe204, re- 
spectively, within the optimized bulk structure and using 
both pure GGA and GGA-f C/. Both the total DOS per 
formula unit and the projected DOS per ion for the d 
states of the various transition metal cations are shown, 
the latter separated into t2g/eg and e/t2 contributions, 
respectively. The DOS calculated within LSDA+?/ (not 
shown) do not show any significant differences compared 
to the ones calculated using GGA+U . 

From the projected DOS it can be seen that all tran- 
sition metal cations are in high spin states, with one 
spin-projection completely occupied, and that the cu- 
bic component of the crystal field on the octahedrally- 
coordinated [Oh) sites lowers the t2g states relative to 
the Cg states, whereas on the tetrahedrally-coordinated 
[Td) sites the e states are slightly lower in energy than the 
t2 states. Due to the Neel-type ferrimagnetic order, the 
local majority spin direction on the sites is reversed 



FIG. 3: (Color online) Total and projected DOS per formula 
unit for CoFe2 04. Left (right) panels correspond to GGA 
(GGA-)~(7) calculations. The d states of Co{Oh) (upper pan- 
els), Fe(Oh) (middle panels), and Fe(Td) (lower panels) are 
separated into t2g (green/dark grey) and Cg (red/black) con- 
tributions for the Oh sites and into e (green/dark grey) and t2 
(red/black) contributions for the Td sites. The total DOS is 
shown as shaded grey area in all panels. Majority (minority) 
spin projections correspond to positive (negative) values. 
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relative to the Oh sites. 

It is apparent that within GGA CoFe204 turns out 
to be a half-metal, in contrast to the insulating behav- 
ior found in experiment. ^"^ This is similar to what has 
been found in previous LSDA calculations.^ The half- 
metallicity is due to the partial filling of the minority 
t2g states of Co(Oh), which in turn results from the for- 
mal configuration of the Co^+ cation (see upper left 
panel of Fig. [S]). This apparent deficiency of the GGA 
approach is corrected within the GGK+U calculation as 
can be seen in the right part of Fig. |3l We note that 
this is very similar to the case of rocksalt CoO, which 
also contains Co^"*" with a dJ electron configuration that 
leads to a metallic solution in pure LSDA,*^ whereas ap- 
plication of the DFT-f C/ approach leads to an insulating 
state in agreement with the experimental observations/^ 



FIG. 4: (Color online) Total and projected DOS per formula 
unit for NiFe2 04. Left (right) panels correspond to GGA 
(GGA+C/) calculations. The d states of Co(Oh) (upper pan- 
els), Fe(Oh) (middle panels), and Fe{Td) (lower panels) are 
separated into t2g (green/dark grey) and Cg (red/black) con- 
tributions for the Oh sites and into e (green/dark grey) and t2 
(red/black) contributions for the Td sites. The total DOS is 
shown as shaded grey area in all panels. Majority (minority) 
spin projections correspond to positive (negative) values. 



TABLE IIL Calculated magnetic moments (in ^b) for bulk 
CoFe2 04 and NiFe2 04 using different exchange-correlation 
functionals. The total magnetic moment per formula unit 
amounts to 3 (2 hb) for CoFe204 (NiFe204), respectively. 
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E-Epei-mi [eV] 



6-4-2 2 



The Hubbard correction splits the occupied and unoc- 
cupied parts of the minority spin t2g states of the Co^+ 
cations (see upper right panel of Fig. [3]) , thereby opening 
an energy gap. In addition, the local spin splitting on 
the Fe cation is drastically enhanced, shifting the local 
majority spin d states towards the bottom of the valence 
band. We note that once the value of Ucs on the Co sites 
is large enough to push the corresponding unoccupied 
minority spin t2g states above the lowest minority-spin 
Fe{Oh) states, the width of the band gap is determined 
by the difference in energy between these lowest unoc- 
cupied minority-spin Fe states and the highest occupied 
minority spin t2g states of Co(0/i). Therefore, a fur- 
ther increase of UcS on the Co sites does not significantly 
change the size of the band gap. Similarly, the band gap 
depends only weakly on the specific value of UcS on the 



CoFe2 04 


Co(0;,) 


Fe(Ofe) 


Fe(Td) 


LSDA-HU 


+2.52 


+3.99 


-3.82 


GGA 


+2A3 


+3.66 


-3.45 


GGA-HU 


+2.62 


+4.10 


-3.98 


NiFe2 04 


Ni(Oft) 


Fe(O,0 


Fe(Td) 


LSDA-f-U 


+1.49 


+4.00 


-3.82 


GGA 


+1.36 


+3.71 


-3.46 


GGA-fU 


+1.58 


+4.11 


-3.97 



Fe sites. 

The gap size of 0.9 eV obtained within GGA+C/ for 
the chosen values of UcS is comparable to the 0.63 eV ob- 
tained by Antonov et al. using LSDA+C/ with C/off=4.0 
eV for the Co(0,,) and C/cff=4.5 eV for the Fe{Oh) and 
Fe(Trf) cationS)^ and also agrees well with the value 
of 0.8 eV reported by Szotek et al. utilizing a self- 
interaction corrected LSDA approachi^ 

Ni^+ formally has one additional electron compared to 
Co^+, leading to a fully occupied minority spin t2g mani- 
fold for the Ni^"*" cation within a cubic (Oh) crystal field. 
Accordingly, NiFe204 exhibits a tiny gap of '^0.1 eV be- 
tween the occupied minority spin t2g states of Ni(0;i) 
and the unoccupied minority t2g states of Fe(0;i) even 
in pure GGA (see left panels of Fig. |4]). However, the 
use of GGA+f/ leads to a significant enlargement of this 
energy gap to a more realistic value of 0.97 eV for the 
chosen values of Ucff. This is in good agreement with 
band gaps of 0.99 eV and 0.98 eV reported by Antonov 
et al^^ and Szotek et al.,^^ respectively. Similar to the 
case of CoFe204 the Hubbard correction also leads to a 
strong enhancement of the local spin splitting on the Fe 
sites in NiFe2 04. 

Table IIIII shows the local magnetic moments of the 
transition metal cations per formula unit calculated 
within the three different approaches. The total mag- 
netic moment is independent of the applied exchange- 
correlation functional, and equal to the integer value that 
follows from the formal electron configuration of the tran- 
sition metal cations and the Neel-type ferrimagnetic ar- 
rangement (3 for CoFe204 and 2 /is for NiFe2 04). 
This is a result of the either half-metallic or insulating 
character of the underlying electronic structures. Never- 
theless, as can be seen in Table Hill the Hubbard correc- 
tion has a significant infiuence on the spatial distribution 
of the magnetization density and the use of GGA+C/ (and 
LSDA+C/) leads to more localized magnetic moments 
compared to GGA, indicated by the increased magnetic 
moments corresponding to the different cation sites. 

The results presented in this section indicate that for 
a realistic and consistent description of the structural, 
electronic, and magnetic properties of both CoFe204 and 
NiFe2 04, a Hubbard correction to either LSDA or GGA 
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FIG. 5: (Color online) Calculated c/a ratio of NiFe204 as a 
function of epitaxial strain Exx obtained from GGA+C/ calcu- 
lations and different cation arrangements ("settings") on the 
Oh sites. 



TABLE IV: Bulk modulus B, 2-dimensional Poisson ratio 
1^20, and elastic coefficients Cn and C12 for CoFe2 04 and 
NiFe2 04, obtained for different exchange-correlation poten- 
tials and different cation arrangements ("setting" s), in com- 
parison to experimental data. The experimental 1/20 has been 
evaluated from Eq. ((7| using the experimental elastic con- 
stants. 



CoFe2 04 
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B [GPa] 


1^20 


Cii [GPa] 


C12 [GPa] 


LSDA+U 


1 


206.0 


1.191 


282.0 


167.9 




3 




1.185 


282.7 


167.6 


GGA+U 


1 


172.3 


1.147 


240.8 


138.1 




3 




1.132 


242.5 


137.3 


Exp.i^ 




185.7 


1.167 


257.1 


150.0 


NiFe2 04 


S 


B [GPa] 


'^2D 


Cii [GPa] 


C12 [GPa] 


LSDA+C/ 


1 


213.1 


1.172 


294.4 


172.5 




3 




1.167 


295.1 


172.2 


GGA+U 


1 


177.1 


1.115 


251.2 


140.0 




3 




1.106 


252.2 


139.5 


Exp.*2 




198.2 


1.177 


273.1 


160.7 



is required. In the following we will therefore present 
only results obtained within the LSDA+[/ and GGA+U 
approaches. 



B. Epitaxial strain and elastic properties 

Fig. [5] shows the relaxed c/a ratio of NiFe2 04 as 
function of the epitaxial constraint Sxx, obtained from 
GGA+U calculations as described in Sec. Ill CI The case 
of CoFe204 is very similar. Two important things can 
be seen from this. First, due to the orthorhombic Imma 
symmetry of the chosen cation arrangement on the Oh 
sites, the c/a ratio is not exactly equal to 1 at zero strain. 
However, this effect is clearly negligible compared to the 
changes in c/a induced by epitaxial strains of order '^1 % 
and therefore docs not affect our further analysis. Sec- 
ond, the slope of the c/a ratio, which characterizes the 
elastic response of the material, is nearly completely un- 
affected by the different cation arrangements. 

From the data shown in Fig. [5] we can therefore ob- 
tain the 2-dimensional Poisson ratio 1^20 (Eq. ([7])), which 
relates in-plane and out-of-plane strains. Together with 
the bulk moduli listed in Table HI we can then determine 
the two elastic constants Cu and C12 from Eqs. ^ and 
©• 

The calculated 2-dimensional Poisson ratios 1/20 and 
elastic constants Cu and C12, together with the bulk 
moduli already presented in Table HI are listed in Ta- 
ble IIVI and are compared with experimental results from 
Ref. 42. It can be seen that, as already pointed out, the 
specific cation arrangement has nearly no influence on the 
value of 1^20 and thus Cu and C12. On the other hand, 
the specific choice of either LSDA-|-?7 or GGA+U has a 
noticeable effect. Similar to the case of the bulk modulus, 
we observe an overestimation (underestimation) of the 
elastic constants Cn and C12 in the LSBA+U {GGA+U) 



FIG. 6: (Color onhne) CoFe204 GGA+U total energy dif- 
ferences for orientation of the magnetization along various 
crystallographic in-plane directions with respect to the [001] 
direction (black solid lines): A [100], T [010], M [110], and 
► [110]. Red broken lines denote crystallographic directions 
which include also out-of-plane components, namely [101] 
(dashed-dotted line), [Oil] (dashed-double dotted line), and 
[111] (dashed line). Left (right) panels contain the results 
corresponding to setting 1 (setting 3). 
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calculations. The same general trend holds for the 2- 
dimensional Poisson-ratio of CoFe204, while for NiFe204 
the use of LSDA-I-C/ also slightly underestimates 1^20- 
Overall the deviations are only within a few percent of 
the experimental data (1-2 % for CoFe204 and up to 
6 % for NiFe2 04), and we therefore conclude that both 
LSDA-|-[/ and GGA-I-J7 allow for a good description of 
the strain response of CoFe204 and NiFe2 04. 



C. Magnetoelastic coupling 

The calculated MAEs, defined as the energy differ- 
ences for various orientations of the magnetization rel- 
ative to the energy for orientation of the magnetization 
parallel to the [001] direction, are depicted in Fig. [5] for 
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FIG. 7: (Color online) NiFe204 LSDA+(7 (GGA+17) total 
energy differences for orientation of the magnetization along 
various crystallographic in-plane directions with respect to 
the [001] direction are shown in the upper (lower) panels 
(black solid lines): A [100], T [010], < [110], and ► [110]. Red 
broken lines denote crystallographic directions which include 
also out-of-plane components, namely [101] (dashed-dotted 
line), [Oil] (dashed-double dotted line), and [111] (dashed 
line). Left (right) panels contain the results corresponding to 
setting 1 (setting 3). 
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CoFe204 (GGA+[/ only) and m Fig.[7]for NiFe204 (both 
LSDA+f/ and GGA+C/). It can be seen that the calcu- 
lated MAEs for CoFe204 are roughly five to six times 
larger than in NiFe204. Furthermore, to a good approxi- 
mation, the calculated energy differences exhibit a linear 
dependence on strain. Deviations from this linear behav- 
ior are most pronounced for the case of NiFe2 04 within 
setting 3. Since the pure elastic response shown in Fig. El 
does not exhibit any significant non-linearities, we con- 
clude that higher order magnetoelastic terms are respon- 
sible for the slightly non-linear behavior of the MAE in 
this case. 

It can also be seen, that in all cases the strain depen- 
dence, i.e. the slope of the various curves shown in Figs.lB] 
and [71 is largest for the in-plane versus out-of-plane en- 
ergy differences, i.e. for orientation of the magnetization 
along the [100], [DIG], [110], and [110] directions (com- 
pared to the [001] direction), consistent with Eqs. ([5]). 
Thus, tensile strain favors perpendicular anisotropy (easy 
axis perpendicular to the "substrate") and compressive 
strain favors in-plane orientation of the magnetization, 
i.e. Bi > 0. For sufficient amount of strain the easy axis 
of magnetization will therefore always be oriented either 
in-plane or out-of-plane, consistent with various experi- 
mental observations in thin GoFe204 films under tensile 
strain^^^ 

According to the phenomenological magnetoelastic 



theory for a cubic crystal discussed in Sec. Ill Bl in partic- 
ular Eqs. the strain dependence of the in-plane versus 
out-of-plane anisotropy should be stronger by a factor of 
2 compared to the anisotropy corresponding to [101] or 
[Oil] orientation of the magnetization, and by a factor of 
3/2 compared to [111] orientation. We note that these 
ratios are very well observed by the calculated anisotropy 
energies shown in Figs. [Hand [71 This indicates that the 
strain-dependence of the calculated anisotropy energies 
is rather independent of the specific cation arrangement, 
which allows us to obtain the magnetoelastic constants 
of CoFe204 and NiFe204 from our calculations. 

On the other hand, we recall that due to the spe- 
cific cation arrangement used in the calculations, and the 
resulting symmetry lowering from cubic to orthorhom- 
bic, the three cubic axes are not equivalent even for 
zero strain. For setting 1 the y and z directions are 
equivalent, but different from a;, whereas for setting 3 
X and y are equivalent but different from z. This is 
reflected in the calculated anisotropy energies for zero 
strain (exx — 0), which in setting 1 are largest between 
z/y and the x direction. (Note that for zero strain both 
settings are completely equivalent apart from a rotation 
of the coordinate axes by 120° around the [111] direc- 
tion.) The anisotropy induced by the symmetry-lowering 
to orthorhombic is therefore at least of the same magni- 
tude as the cubic anisotropy in the disordered inverse 
spinel structure, and the calculated anisotropy energies 
for zero strain are therefore not representative for the in- 
verse spinel structure with random distribution of cations 
on the Oh sites, i.e. the cubic anisotropy constant K can- 
not be determined from our calculations. 

From the data shown on Figs. [Hand [71 we thus obtain 
the magnetoelastic coefficients Bi by using the appropri- 
ate equation out of Eqs. ^ (together with the calculated 
values for ^"20) for each calculated energy difference in- 
dividually, and then average over the resulting values for 
Bl within the same setting and for the same exchange- 
correlation functional. From the so-obtained magnetoe- 
lastic coefficients Bi we then calculate the linear magne- 
tostriction coefficient Aioo from Eq. ((H]) using the elastic 
moduli determined in Sec. IIII Bl The resulting values for 
both Bl and Aioo are listed in Table IVl 

It can be seen that overall the calculated magnetostric- 
tion constants are in very good agreement with avail- 
able experimental data, despite the difficulties related to 
the specific choice of cation arrangement and exchange- 
correlation functional. In particular the large difference 
in magnetostriction between CoFe204 and NiFe2 04 is 
well reproduced by the calculations. For NiFe204 the 
difference in the calculated values for Aioo due to the 
use of either LSDA-I-J7 and GGA+U is larger than the 
effect of the different cation settings. For CoFe204 the 
cation arrangement seems to have a larger influence on 
Aioo than for NiFe2 04. This is consistent with the large 
spread in the experimentally obtained magnetostriction 
for CoFe204, where different preparation conditions can 
lead to differences in cation distribution/inversion or 
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TABLE V: Magnetoelastic coupling constant Bi and mag- 
netostriction constant Aioo for CoFe204 and NiFe204 using 
different exchange-correlation functionals and cation arrange- 
ments ("setting" s) in comparison with available experimental 
data. 
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(xlO^'*) 


LSDA+U 


1 






10.03 


-54.9 




3 






9.65 


-52.3 


GGA+U 


1 


29.23 


-189.7 


6.72 


-40.3 




3 


39.74 


-251.7 


6.08 


-35.9 


Exp. 






-225.0" 
-250.0' 
-590.0" 




-50.9" 
-36.0" 



"Polycrystallinc CoFe204 (Ref.,49). 

'Single crystals with C01.1Fe1.9O4 composition (Ref. . 50^ . 
"Single crystals with Coo.8Fe2.2O4 composition (Ref. Isd). 
"^Single crystals of NiFe2 04 (Ref.^l). 

"Single crystals with Nio.8Fe2.2O4 composition (Ref. [soh . 



slightly ofF-stoichiometric compositions. In addition, only 
very few measurements have been performed on single 
crystals, whereas for polycrystalline samples only a su- 
perposition of Aioo and Am is measured. 

The rather good agreement between our calculated val- 
ues for Aioo and the available experimental data demon- 
strates that in principle a quantitative calculation of 
magnetostriction in spinel ferrites is feasible, in spite of 
the difficulties related to the inverse spinel structure with 
random cation distribution on the B site, and the usual 
difficulties regarding an accurate description of exchange 
and correlation effects in transition metal oxides. 



IV. SUMMARY AND OUTLOOK 

In summary, we have presented a systematic first prin- 
ciples study of the effect of epitaxial strain on the struc- 
tural and magneto-structural properties of CoFe204 and 
NiFe2 04 spinel ferrites. Special care was taken to as- 
sess the quantitative uncertainties resulting from differ- 
ent treatments of exchange-correlation effects and dif- 
ferent cation arrangements used to represent the inverse 
spinel structure. 

It has been shown, in agreement with earlier works, 
that "beyond LSDA/GGA" methods are required for a 
proper description of the electronic and magnetic proper- 
ties of spinel ferrites CoFe2 04 and NiFc2 04. The "+f7" 



approach used in the present work leads to a realistic elec- 
tronic structure and good quantitative agreement with 
available experimental data for lattice, elastic, and mag- 
netostrictive constants. 

We find that the specific cation arrangement used to 
represent the inverse spinel structure has only little effect 
on the structural properties. The corresponding effect on 
the magnetoelastic constants is also weak, with a some- 
what stronger influence in the case of CoFe2 04. The 
latter fact is consistent with the considerable spread in 
the reported values for the magnetostriction constant for 
different samples of this material. In general, the starting 
point for the "+U" correction, i.e. either LSDA or GGA, 
has a somewhat stronger influence on the calculated ma- 
terials constants than the different cation settings. 

However, in spite of these uncertainties, the overall 
agreement between our results and experimental data is 
very good. In particular, the calculated magnetostric- 
tion constants Aioo of about —45 ppm for NiFe2 04 and 
about —220 ppm for CoFe204 fall well into the spectra 
of available experimental values obtained from different 
samples (see Table IVT) . Consistent with the negative sign 
of Aioo, the easy magnetization direction changes from 
in-plane for compressive epitaxial strain to out-of-plane 
for tensile strain. This gives further confirmation that 
the reorientation of the easy axis observed experimen- 
tally in thin films of CoFe204 under different conditions 
is indeed predominantly strain-drivenJ^ii^i^ 

In summary, our results indicate that a quantitative 
description of both structural and magnetoelastic prop- 
erties in spinel ferrites is possible within the DFT+U 
approach, which opens the way for future computational 
studies of these materials. Such calculations can then 
provide important information regarding the effect of 
cation inversion and off-stoichiometry, which can be used 
to optimize magnetostriction constants and anisotropy in 
spinel ferrites. 
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